Segmented thermoelectric generators: Impact of junction temperature variation on 

the series circuit's properties 



(N 

o 

G 

in 



CO 



C 

O 

o 



> 

OS 

in 

O 

(N 



X 



Y. Apertet, 1 ^ H. Ouerdane, 2 C. Goupil, 3 and Ph. Lecoeur 1 

1 Institut d'Electronique Fondamentale, Universite Paris-Sud, CNRS, UMR 8622, F-91405 Orsay, France 
2 CNRT Materiaux UMS CNRS 3318, 6 Boulevard Marechal Juin, F-14050 Caen Cedex, France 
3 Laboratoire CRISMAT, UMR 6508 CNRS, ENSICAEN et Universite de Caen Basse Normandie, 
6 Boulevard Marechal Juin, F-14050 Caen, France 
(Dated: March 8, 2013) 

We present and analyze a model of a system constituted of two thermoelectric generators (TEGs) 
thermally and electrically connected in series. For simplicity the parameters of both TEGs are 
assumed to be temperature-independent. The derivation of the equivalent electrical series resistance 
yields an unexpected term whose meaning and impact are discussed. We demonstrate that this term 
must exist as a consequence of thermal continuity at the interface, since it is related to the variation 
of the junction temperature between the two TEGs as the electrical current varies. We then derive an 
expression for the equivalent series figure of merit and use it to discuss the approach of thermoelectric 
compatibility. Finally we highlight the strong thermal/electrical symmetry between the parallel and 
series configurations. 

PACS numbers: 84.60.Rb, 72.20.Pa, 84.60.Bk 



I. INTRODUCTION 

Optimization of thermoelectric systems for energy con- 
version applications not only involves the improvement of 
the materials' properties to enhance the so-called figure 
of merit, but also a strategic reflexion on device design 
[l| . The segmentation of thermoelectric legs in a module 
is one of the strategies used to yield the best possible per- 
formance of the materials. A segmented leg is obtained 
by stacking different thermoelectric materials in such a 
way that each part of the leg is optimized for the tem- 
perature that it experiences within the structure. The 
basic principle of this technique lies in the temperature- 
dependence of the materials' properties: a given material 
presents interesting conversion capabilities only within a 
limited temperature range. So, when the figure of merit 
of a given material collapses at some point of the leg be- 
cause of the variation of temperature along the leg, this 
material is replaced at this point by another one, bet- 
ter suited to the local temperature. The benefits on the 
global performance of the device is obvious when the tem- 
perature difference imposed on the TEG is important Q ■ 
A high performance leg is composed of different segments 
thermally and electrically in series. 

The principle of segmentation was patented in 1962 
However this segmentation strategy is not always 
the best solution to achieve high power efficiency as it 
was pointed out many times and even before the publi- 
catioon of the patent. Thus, along the years, many cri- 
teria were selected and used to see if the global efficiency 
could sufficiently increase by segmentation. If some stud- 
ies focused on local entropy production [J] and differen- 
tial output power [j| , others developped more global ap- 
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FIG. 1. Schematic representation of a segmented thermoelec- 
tric generator composed of two parts. 



proaches @, 0] • Recently Snyder and Ursell proposed the 
so-called compatibility approach giving a rule of thumb 
based on a compatibility factor associated with each ma- 
terial [8|, [9j . The associated thermodynamic framework 
was further discussed in Ref. 10'. Our purpose here is to 
highlight the underlying physical mechanisms permitting 
or not compatibility between different thermoelectric ma- 
terials. We focuse on the simple case of a leg composed 
of only two thermoelectric segments. Even if the method 
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is intended to be applied to materials with temperature- 
dependent parameters, for the sake of simplicity in this 
article, we develop and analyze a model with constant 
parameters. 

Our paper is organized as follows. In Sec. [HI we present 
our model of a two-sgement thermoelectric generator; we 
focus especially on the equivalent series parameters and 
figure of merit and on their physical meaning. Then in 
Sec. IIIII we compare our derivations to the so-called com- 
patibility approach of the segmented generators, demon- 
strating that the latter works accurately only under re- 
strictive conditions. In Sec. IIV1 we relate the system 
treated in this article to a device composed of two TEGs 
thermally and electrically in parallel: we highlight the 
fact that there is a surprising thermal/electrical symme- 
try between the two configurations. We end this paper 
with concluding remarks and an appendix where we show 
that the neglect of Joule heating and transferred power 
to the load in our model, is of no importance in the cal- 
culation of the junction temperature. 



II. EQUIVALENT MODEL OF TWO 
THERMOELECTRIC GENERATORS IN SERIES 



■Kconv.i and i^ C onv,2- These conductances depend on elec- 
trical load through the electrical current I and may be 
expressed as [L?T ]: 

a{TIi 

^conv,i — (o) 

This model is derived assuming that the heat flux is con- 
stant along a TEG and equal to its average value. This 
assumption amounts to neglect both the Joule heating 
and the electrical power produced, i.e., considering a 
TEG with small efficiency. Here, we further make the 
approximation that the mean temperature T can be con- 
sidered as constant through the system as it is high com- 
pared to the temperature difference between the heat 
reservoirs. All these assumptions are consistent with the 
linear approximation used to describe the behavior of the 
segments. 



B. Temperature at the junction 

Continuity of heat flux at the junction yields: 



A. Definitions 

We consider two thermoelectric generators labeled 1 
and 2, electrically and thermally connected in series as 
depicted in Fig. [1] The electrical current through the 
resistive load i?i oa d, is /. Each generator is character- 
ized by an internal electrical resistance Ri, a Seebeck 
coefficient on, and a thermal conductance under open- 
circuit condition Ki (i=l,2). The voltages and tempera- 
ture differences across each module are denoted AV* and 
AT;. The electrical currents and average thermal fluxes 
through each TEG, U and Iq 4 , are derived from the phe- 
nomenological force- flux formalism fTTj | : 
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(1) 



where T is the average temperature of the whole system. 
The figure of merit of each TEG is Z{T = c$T/{RiKi). 

The key issue that we address in this paper is the deter- 
mination of the temperature at the junction, T m . This 
topic was hardly considered even though El-Genk and 
Saber noticed that the electrical current could have an 
impact on this temperature [2j. 

The mean heat flux inside a segment, i.e., in one of the 
TEGs, is given by: 



KkCThot - T m ) + ai TI = K 2 (T m - T cold ) + a 2 TI (4) 

from which we derive the expression for the temperature 
at the junction between the two segments, T m : 
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K x + K x 



[(ax - a 2 )TI + if iT hot + K 2 T coM ] (5) 



For a more accurate description of the heat flux conti- 
nuity, the reader may refer to the appendix, where we 
show that this simple approach gives results pretty close 
to those obtained without approximation. 

Let us now focus on the eletrical part of the generator. 
The voltage across the whole module is the sum of the 
voltages across each segment: 



AV = ai (T hot - T m ) + a 2 (T m - T coM ) -(Rx+R 2 )I (6) 
Substitution of T m by its full expression (|5|) yields 



K 2 on, + Kt. ceo , 
AV = V „ 2 (T hot - T cold ) 



Kx+K 2 
{ ai ~a 2 ) 2 T 
K x + K 2 



Rx + R 2 



(7) 



ctiTIi + KiATi 



(2) 



In Fig. [T] the heat conveyed by the electrical current 
is represented by the additional thermal conductances 



and we thus obtain an expression with the same form as 
the one expected for a Thevenin generator: 



AV = a eq (Thot - T cold ) - R cq I 



(8) 



where a oq is the equivalent series Seebeck coefficient and 
i? eq is the equivalent series electrical resistance. By iden- g~ 
tification we get the following explicit expressions for 
both equivalent quantities: 



K 2 ai + K\ol-x 
K x + K 2 



R eq = Rl + R2 + Rr, 



(9) 



(10) 



The appearance of a finite third term in the equation 
above is rather unexpected but, as we will explain below, 
it is a crucial term; its expression reads: 
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When considering the system in the open circuit consid- 
eration we find that the equivalent thermal conductance 
K eq for the whole system is given by 



K K 2 K~i 
cq ~ K 1+ K 2 > 



(12) 



which is a standard form for an equivalent series conduc- 
tance. 



C. On the meaning of i? re iax 

The quantity -Rreiax is proportional to the difference be- 
tween the Seebeck coefficients a\ and a 2 , which clearly 
shows that in presence of materials with significantly dif- 
ferent thermopowers, the standard form for the equiva- 
lent series resistance does not apply; further, note that 
^reiax is by no means related to an electrical contact re- 
sistance at the interface between the two segments, which 
we neglect in our model. Its appearence in Eq. (|10[) is 
directly related to the relaxation of the temperature at 
the junction T m , which depends on the value of the elec- 
trical current I. Indeed, to ensure thermal continuity at 
the interface, when the convective part of the thermal 
current varies, the conductive part must change accord- 
ingly to account for the modification of the temperature 
at the junction. 

Thus when for exemple ot\ < a 2 , as the electrical cur- 
rent is constant along the whole system, the conductive 
part of the thermal flux is larger in the second segment 
than it is in the first one. Thus, to satisfy the condition of 
continuity of the thermal flux, the conductive part must 
be larger in the first segment than it is in the second 
one. This thermal balance can only be obtained by a 
relaxation of the temperature T m . Since this difference 
between the convective thermal fluxes in each segment 
increases when the electrical current increases, T m varies 
with the electrical current /, as displayed in Eq. ([5]). The 



FIG. 2. Temperature difference at the edges of each segment 
when the global system is at maximum efficiency as a function 
of ai (normalized by 0:2). TEG\ in full line and TEGi in 
dashed line 



result of such temperature relaxation at the junction al- 
ways leads to increase the temperature difference at the 
edges of the segment exibiting the lowest Seebeck coeffi- 
cient, i.e., the one with less capability to produce power, 
and conversely to increase the temperature difference at 
the edges of the segment exhibiting the highest Seebeck 
coefficient: the global electromotive force developped by 
the whole system is thus reduced. This detrimental ef- 
fect for energy conversion performances is reflected by 
the existence of the additional term -R ro iax- 

To further discuss the consequences of the presence 
of -R re iax in Eq. (fT0|) . we propose a numerical example. 
The parameters for TEG 2 are set to the following val- 
ues: Z 2 T = 1, K 2 = 2.5 mW/K, R 2 = 4.8 mfi, and 
a 2 = 200 /iV/K. These parameters are consistent with 
the properties at room temperature of a particular bis- 
muth telluride compound, (Bio.25Sbo.75 VTe3, one of the 
most efficient thermoelectric materials [1J|. For TEG 1, 
we make the assumption that the material used has a 
fixed figure of merit Z\T = 2, which is optimistic but 
remains feasible [3]. K\ is also fixed, with K\ = K 2 . 
The thermopower a± is variable and so is R\ since it is 
matched to the value of Z\T . The reservoirs tempera- 
tures are Thot = 305 K and T co id = 295 K, and the mean 
temperature is T = 300 K. 

Figure © displays the dependence on a± of the tem- 
perature difference experienced by each module when the 
whole system works at maximum efficiency. The thermal 
conductances of each TEG are identical; this specifica- 
tion allows to obtain an equal partition of the tempera- 
ture difference when there is no electrical current inside 
the structure. When the Seebeck coefficients are identi- 
cal there is, as expected, no relaxation of T m , but when 
ot\ differs from a%, we notice that one segment experi- 
ences a larger temperature difference. The favored side 
is always the one with the smaller Seebeck coefficient. 



D. On the equivalent series figure of merit 

Now that we have derived the equivalent parameters 
we may express the equivalent series figure of merit Z eq 
as follows: 



cq 



Zcq ~ R K 



(13) 



which, using Eqs. (|9|). ()10|) and fj 12[) . may be rewritten as 
the product of two terms: 



Z^ q — Y -^series 



(14) 
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FIG. 3. Z cq T, Z SCI i cs T and Y plotted against the ratio 01/02- 
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and 



Y 



(ai - a 2 ) T 
{R l +R 2 ){K 1 +K 2 ) 



(16) 



The term Series should be seen as an expected term for 
which the equivalent series electrical resistance is given 
by the standard sum of R\ and R 2 . The factor Y, on 
the contrary, is purely related to the relaxation of the 
temperature at the junction between the two TEGs. This 
factor reflects the fact that this relaxation always leads to 
a decrease of the performance; it is indeed smaller than 
1 except for ct\ = a 2 . It is also interesting to note that 
the thermal conductances and electrical resistances also 
appear in the expression of Y and that these quantities 
must be small for Y to have a significant impact on Z eq . 

Figure [3] obtained for the same numerical parameters 
as above, illustrates the fact that in order to optimize 
Z cq , one has to make a compromise between the mere 
association of the segments, represented by Z ser i eSl and 
the effect resulting from the mismatch between the See- 
beck coefficients, represented by Y. The attenuation of 
the performances due to the thcrmopower mismatch calls 
to mind the concept of thermoelectric compatibility dis- 
cussed by Snyder and Ursell [9(]. 

To end this section on the equivalent series figure of 
merit, we want to stress that Bergman and Levy's the- 
orem is satisfied: Z cq T, given by Eq. (fT4")) . is always 
smaller than that of the segment with the highest fig- 
ure of merit [HI]. In the next section we compare our 
treatement of segmentation to that of the compatibility 
approach. 



III. ON THE THERMOELECTRIC 
COMPATIBILITY 

The optimization of segmented TEGs is inseparable of 
the concept of thermoelectric compatibility developped 
by Snyder and Ursell @. We show in this section that 
our approach is, to some extent, indeed different from the 
compatibility approach, and we explain why we observe a 
discrepancy between the results given by both methods. 

To evaluate the difference between the two approaches, 
we first remind the definitions of the physical quantities 
pertaining to the thermoelectric compatibility. The rela- 
tive current it, for TEG i is defined as: 



(17) 



where J is the electrical current density and m the 
thermal conductivity under open-circuit condition. The 
macroscopic expression for the relative current is: 



Ui 



KiATi 



(18) 



The compatibility factor s, namely the optimal value 
for the relative current u is expressed as: 



y/1 + ZjT - 1 
ctiT 



(19) 



Two materials are considered compatible only if their 
compatibility factors are close (within a factor of 2) [8J. 
This condition can be viewed as a necessity to have op- 
timal values of the electrical currents in each segment, 
which are close enough: in this way, as the current is the 
same through the device, each segment works optimally 
at the same time (l6| . 

Our method is based on the optimization of the equiv- 
alent series figure of merit Z eq T defined in Eq. (Ti"4"|) and 
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satisfaction of the condition of heat continuity (though, 
still neglecting the produced power and the Joule effect 
at the interface) the conductive part given by KiATi is 
internally set to the appropriate value through the re- 
laxation of T m . The compatibility approach thus looses 
accuracy when the figures of merit of each segment are 
both high (i.e. greater than 1) and relatively dissimi- 
lar, which is the case of the numerical example treated 
in this paper. The determination of the relative current 
at a local scale through numerical simulations remains a 
powerful tool for device optimization nonetheless, as it 
allows to deal with generators composed of more than 
two segments and with materials characterized by non 
constant parameters [H, [Til H3] ■ 



FIG. 4. Maximum efficiency (scaled the Carnot efficiency) 
as a function of ot\ evaluated both numerically and from 
Eq. (J20J). 



IV. COMPARISON BETWEEN SERIES AND 
PARALLEL CONFIGURATIONS 



displayed in Fig. |3] for the particular case considered in 
this paper. To check that Z cq T is a meaningful quan- 
tity we evaluate the maximum efficiency using the tradi- 
tional expression containing the figure of merit (see, e.g., 
Ref. [13) with the value of Z cq T: 



^Carnot 



ZcqT 



V 1 + Z ^ T + jfif 



(20) 



The result is displayed in Fig. For comparison 

we also plot the maximum efficiency calculated with a 
numerical simulation using the complete description of 
the device given in Eq. (fTJ) . The two curves are very close 
so we may conclude that Z eq T is a proper parameter to 
estimate the efficiency capability of the whole device. 

With the case considered in this paper, the value of a\ 
that permits satisfaction of the condition s\ — S2, is given 
by a\ = a 2 -^==57 = 1-77 cti (vertical dotted-line in 
Fig. (|4}) whereas with our equivalent model we find an 
optimal value of a\ = 1.34 ai-, which, when compared to 
the numerical result, leads to a better optimization for 
the maximum efficiency. 

Why the compatibility approach fails to accurately 
predict the optimal set of parameters in this particular 
case? To derive their expression of the compatibility fac- 
tor s, it appears that Snyder and Ursell assumed that 
heat transport by conduction remains constant along the 
device: "Since all segments in a thermoelectric element 
are electrically and thermally in series, the same current 
I and similar conduction heat Ak\7T flow through each 
segmen t" |9| . This assumption is more explicitely stated 
in Ref. [16j . However, as soon as the figure of merit varies 
from one segment to the next one, the partition of heat 
flux between the convective and the conductive contribu- 
tions changes along the device. As derived in Ref. fl8j . 
ZT may indeed be related to the ratio of these two con- 
tributions and it follows that, in order to guarantee the 



In a recent article [2l| we studied the association of 
two TEGs electrically and thermally connected in paral- 
lel: we demonstrated that if the Seebeck coefficients of 
the two segments composing the device are different, an 
internal electrical current may develop and yield an ad- 
ditional term for the equivalent parallel conductance un- 
der open circuit condition. This additional term, K com , 
results from the convection process, i.e. heat that is con- 
veyed by the electrical current inside the structure. Note 
that in Ref. [2l[ we used the word advection but we real- 
ized since then that convection is better suited to char- 
acterize this heat transport process. The conductance 



(ai 



a 2 ) 2 T 



i?l + i?2 



(21) 



At this stage, it is instructive to highlight the similar 
forms of Eq. (|2"Tj) and Eq. (jXTJ) : both expressions are pro- 
portional to the square of the difference of the Seebeck co- 
efficients, and one may switch from one to the other only 
by exchanging the nature of the physical property under 
consideration: one only needs to replace the electrical re- 
sistance by the thermal conductance and vice versa. It is 
interesting to point out that the underlying symmetry be- 
tween thermal and electrical transport in thermoelectric 
phenomena is reflected by the correspondence between 
these two equations. Indeed we notice that the electrical 
transport impacts the thermal transport in the parallel 
configuration as much as the thermal transport impacts 
the response of the electrical circuit in the series config- 
uration. 

We pursue the comparison between both cases: as for 
the series configuration, the equivalent parallel figure of 
merit may be expressed as the product of an expected 
term and the factor Y as defined in Section HT1 
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Z n T = Y 

eq 



(Ci + g 2 )( Gi g;+gf 

Ki + K 2 



(22) 



which yields the following definition of the equivalent par- 
allel thermopower and conductance: 



II G\a\ + G2OL2 || , , 

a cq = G +G , 2 — and G\ = d + G 2 (23) 

where for convenience, in the parallel case, we use the 
electrical conductance G(= 1/R) instead of the electrical 
resistance. 



Appendix A: Approximate expression against exact 
expression 

In this article, for simplicity we used an approximate 
description of the TEGs, neglecting both Joule heating 
and power transferred to the load, in order to derive the 
expression of the temperature T m . We show here that 
the approximations we made is justified. To do so, we 
compare the results obtained with Eq. (J5J) and from the 
following exact expression derived from the full condition 
of continuity of the heat flux at the junction: 



l-thot 



KoT, 



2 -L cold 



Ki + K 2 + (a<2 - ai)I 



(Al) 



SUMMARY 



In this article we have presented a simple model of 
a segmented thermoelectric generator composed of only 
two segments. We demonstrated that when the Seebeck 
coefficients of each segment are different, the temper- 
ature at the junction changes as the electrical current 
varies. This effect is embodied in an additional term 
for the equivalent series electrical resistance, and reflects 
the associated decrease of performance. Knowledge of 
equivalent series parameters allowed the derivation of an 
equivalent series figure of merit, which accurately char- 
acterizes the efficiency capability of the whole device. By 
comparison of our approach to the so-called compatibility 
approach we demonstrated that if our simple model sheds 
light on the mechanism of temperature relaxation at the 
junction between each segment, the compatibility factor 
fails to account for this fact. However, at the local scale, 
the relative current remains a suitable and powerful tool 
for efficiency optimization. Indeed it allows consideration 
for devices composed of more than two segments, made 
with materials characterized by temperature-dependent 
parameters. Finally we have highlighted the symmetry 
that exists between the association of TEGs in parallel 
and the association in series. 
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FIG. 5. Comparison of the exact and approximated values of 
T m for the maximum efficiency working condition as a func- 
tion of Qi. 

Retaining the same values of the parameters used for 
the numerical example treated in the main text, figure [5] 
displays the dependence on the thermopower ai of T m , 
given by Eq. (J5J) and by Eq. (jAII) when the whole sys- 
tem works at maximum efficiency. We notice that the 
approximation holds very well for the different values of 
a-y. However, the exact value always is slightly larger 
than the approximated one: this discrepancy is due to 
the Joule heating, neglected in the approximation, which 
raises a little the internal temperature of the structure. 
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